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Abstract 

The problem of target localization with censored noise free binary detectors is considered. In this 
setting only the detecting sensors report their locations to the fusion center. It is proven that if the radius 
of detection is not known to the fusion center, a minimum variance unbiased (MVU) estimator does not 
exist. Also it is shown that when the radius is known the center of mass of the possible target region 
is the MVU estimator. In addition, a sub-optimum estimator is introduced whose performance is close 
to the MVU estimator but is preferred computationally. Furthermore, minimal sufficient statistics have 
been provided, both when the detection radius is known and when it is not. Simulations confirmed that 
the derived MVU estimator outperforms several heuristic location estimators. 
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Target location in the space 
Origin 

Total number of sensors deployed in the region 

Number of detecting sensors 
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Set of locations of detecting sensors 

Vector of locations of detecting sensors 
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A ball with radius R around the i’th sensor 

The possible target region based on Z observation 

The possible target region with sensor jth excluded from evaluation 

Index of sensors forming the minimal sufficient statistic 

Set of locations of the sensors forming the minimal sufficient statistic 

Vector of the locations of sensors forming the minimal sufficient statistic 

Function outputting G\iss from input Z 

Function outputting Z M ss from input Z 

The indicator function of X 

Boundary surface of T 

Probability density function of occurrence of Z if the target located at z T 
Range of random variable X over parameter space 0 
Convex hull of the set Z M ss 

The set of points whose maximum distance from T are not more than R 
A vector storing the distance of (z 1; ..,z k ) elements from z 1 
An estimator based on observation Z 
Center of mass of T(Z) with radius R.\ 
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I. Introduction 

Localization of an unknown transmitter with observations from a network of sensors is a 
well known problem in the literature [HI, (2). The observations can be carried out through 
measurement of Angle of Arrival (AoA) 0], [Q, Time Difference of Arrival (TDoA) [0, 
m, H, or Received Signal Strength (RSS) & ED, DU, GUI, O, d, G3, ED, G2, 
[[H, Ifl9l . [f20j| , [f2TTl . [f22ll . [|23l . When the sensors are mobile, Frequency Difference of Arrival 
(FDoA) can also be used as an additional source of information ll24il . 11251 . Il26l . AoA, TDoA 
and FDoA approaches require sophisticated sensors, and, therefore do not fit well within the 
limitations of wireless sensor networks, specifically for energy and complexity constraints of the 
nodes. In practice, however, an exact(un-quantized) measurement is unrealistic because it requires 
unlimited bandwidth to communicate the data to the fusion center. A binary RSS measurement is 
preferred because it is simpler and requires less resources. A number of papers describe methods 
that use binary RSS localization fl8l . Ifl9l . Il20l . |[2TI . Il22l . Il23l . In some papers it is assumed 
that the propagation model is isotropic and detection is noise free ll20l . [EH, ll22l . If23l . This 
is equivalent to the situation when sensors average the power measurements over a long period 
of time, hence effectively eliminating the effect of noise lf23l . Moreover, the no noise regime 
provides a lower bound for location estimation in the presence of uncertainty such as noise or 
fading. Therefore, a minimum variance unbiased (MVU) estimator of this scenario will serve as 
a benchmark for comparison of the performance of all target localizers. In addition, it provides 
insight into the effect of parameters such as density of sensor deployment and power in the 
localization process regardless of the environment and the type of detection. This assumption 
reduces the detection problem to the question whether the target is located within a certain radius 
of each sensor or not. However, the estimation problem is not well behaved, as the discontinuity 
in the proximity function violates the regularity conditions needed for obtaining Cramer-Rao 
bound (CRB). In addition, some papers consider just the detecting sensors for evaluation of 
the target location lf20l . Il27l . The motivation behind such approach is that when the region of 
sensor deployment is much larger than the detecting radius, the number of non-detecting sensor 
are much higher than the detecting ones. So censoring them will save a lot of communication 
and processing cost ll28ll . while knowing the location of detecting ones still provides a good 
estimate of the target location. 
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In this paper, we first study the redundancy in information provided by detecting sensors and 
find out minimal sufficient statistics for them. Then we investigate the existence of an MVU 
estimator for this problem. Finally, we will introduce some sub-optimal estimators with low 
computational complexity that perform close to optimal and compare the performance of the 
MVU estimator and the sub-optimal estimator with some heuristic ones through simulation. To 
solve the problem in each stage, we have divided the problem into two separate cases: I) when 
the radius of detection is known which is equivalent to the situation when propagation model 
and the transmit power are known to the fusion center; and II) when the radius of detection is 
unknown which is equivalent to the case when propagation model or transmit power are unknown 
such as the case in non-cooperative localization. 

The rest of the paper is organized as follows : section [TT] is the problem formulation, section 
TIT] study sufficient statistics for the observation, section [TV] investigate the existence of MVU 


estimator, section [V] discusses some sub-optimal estimator which behaves close to optimal, 


section VI reports the simulation results and section VII is the conclusion. 


II. Problem Formulation 

Assume that a target is located at unknown location z T = [xi T , x 2t , ..., x Z t ] in l dimensional 
space (in practical applications l is either 2 or 3) and transmits a signal whose power propagates 
isotropically and is attenuated monotonically as a function of distance from the target. N sensors, 
randomly scattered in a deployment region G, with hyper volume A(G). They measure the 
received power and compare it with a threshold, r, to make a binary decision about the target 
presence. We assume the sensors do noise free decision, which can be considered as the limiting 
case when the measured power is averaged over a sufficiently long duration. Furthermore, the 
sensors are configured such that only the detecting sensors report their locations, z 1 ; ..,z n , to 
the fusion center where the localization is performed. Since the received power is a decreasing 
function of distance from the target, there is a ball around the target, Br (z t ), where all the 
sensors inside will detect and those outside will not. From now on we call R the detection 
radius. We assume that G is sufficiently large such that £> 2 ,r(z t ) C G. In addition, we assume 
that at least one sensor detects the target. Let n be the number of detecting sensors (n > 1) and 
S = {1, n} be the set of indices of all detecting sensors. Therefore, Z = {z, : |z e 5} will be 
the set of locations of all detecting sensors and Z = (z.;|i e <S) denote the vector contains those 
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locations. 


III. Sufficient Statistics 

In the this section we derive sufficient statistics for estimation of target location z T . We 
consider the problem in two cases depending on whether the detection radius, R, is known or 
unknown. 

A. Known Detection Radius 

Let us define the possible target region given observation Z as |[29ll . 11301 

T(Z) = n B r ( zO, 

ies 

Alternatively since the order of Z elements does not matter in this definition, we may equivalently 
define T over the Z set i.e. Ii29l . If30l 

T(Z) = f| 8rM, 

ieS 

Also let us define the possible target region with node j removed as 

75 1 z)= n 

ies,i^j 

Clearly, we have z T e T(Z) and T(Z) C Tj(Z ) for all j. Moreover, any sensor whose 
exclusion for evaluating T does not effect the result (i.e. 7}(Z) = T(Z)), does not provide any 
additional information about the target location. Hence, by eliminating all such sensors, we can 
build a new set of indices, 5 M ssQ whose sensors do contribute in shaping the possible target 
region : 

^MSS — oraz)^r(z)}, 

Correspondingly, let Z MSS denote the set of sensors locations whose index are in 5 M ss^ 

^MSS = \ z j\j £ ^mss} , 

and let Z MS s be the vector format of Z MSS i.e. , 

ZmSS = { z j\j £ ^mss) , 

'The subscript MSS is used, since we will proceed to show that .Zmss is in fact the minimum sufficient statistic. 
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We may also represent Z MSS and Z MSS as a function of Z i.e. Z MS s = MSS(Z) or Z MSS = 
MSS(Z). 

Recall that we have assumed that S f 05. Thus, so is <S M ss> and the boundary surface of T, 
denoted by a(T), is composed of hyper spherical domes centered at elements of Z M ss- 

Let 3 be a random variable vector having a probability density function /(Z;z T ) and 0 be 
the z T parameter space. Range of 3 over 0 would be defined as OTll 


ft 0 (3) = {Z|3z T G0; /(Z;0)>O}. 

In this paper, an observation Z is called a possible event if Z e 7?.©(3). 


Lemma 1. For any pair of possible events Z and Z!, 


MSS(Z) = MSS(Z') T(Z) = T(Z'). 


Proof: The forward direction is obvious from the construction of XiSS(Z) and A4SS(Z'). 
For the backward direction let’s assume that it is not true. Then, 


3X,Y e 7Z e (3) , MSS(X) f MSS(Y) & T(X) = T(Y) 


Thus, there is at least one uncommon element between A4SS(X) and A4SS(Y). Without loss 
of generality we assume it to be z p <E MSS(X) and z p f A4SS(Y). Let a represents the 
surface of a set. a(B RjZp ) D cr(T(X)) is a segment of a hyper sphere dome and l + 1 points can 
be selected on it such that they are not located in the same hyper plane. Since T(X) = T(Y), 
<j(Br , Zp ) H a(T(Y)) is also a segment of the same dome and contains those l + 1 points. But 
any l + 1 points on a dome will determine exactly one sensor location belonging to A iSS(Y) 
which has R distance from them all. Thus, z p e A4SS(V) which contradicts our assumption 
and the proof is complete. ■ 

Theorem 1. Z MSS is a minimal sufficient statistic for the estimation of target location. 

Proof: Due to Lemma [lj it suffices to show that T(Z) is a minimal sufficient statistics for 
the estimation of target location. The probability density function that the detecting sensor will 
be located at, Z, can be described as 



z T e T(Z) 
Z T i T(Z) 


( 1 ) 
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where n(Z) represents the number of elements of vector Z and A{.) represents the hyper volume 
of the shape (or simply the area in two dimensional space). Equation ([[]) can be rewritten as 
following which permits the Neyman-Fisher factorization 


= h(Z)g(z T ,T(Z)), 


z T eT(Z) 


where l x is the indicator function of X, h(Z) = 1 ^ and g(z T , T(Z)) 

IztgT(z)- Thus, T(Z) is a sufficient statistics for estimation of z T . 


To prove that it is minimal we consider another sufficient statistics U(.,.) and show that T is 
a function of U through showing that 

V Z ls Z 2 G 7^e(3), W(Z 1 )=W(Z 2 )=^r(Z 1 ) = r(Z 2 ) (2) 

where 3 represents the random variable vector for observations and Z| and Z 2 are two instances 
of that random variable vector. Assume that ([2]) does not hold, then 

3 Zi, Z 2 G 7^e(3); T(Zi) ^ T(Z 2 ) and U(fLf) = U(Z 2 ).eq : minimalU (3) 

Assume that z T G ((T(Zi) — T(Z 2 )) U (T(Z 2 ) — T(Zi))). Without loss of generality we assume 
that z T G T(Zf) — T(Z 2 ). Then /(Z^zx) 0 and /(Z 2 ;z T ) = 0. 

However, due to Fisher Neyman factorization we have 

/(z i;ZT ) = /1«(Z 1 )^( ZT ,W(Z 1 ))^0 

/(Z 2 ;z T ) = /iw(Z 2 ) 5 ^(z T ,W(Z 2 )) = /i«(Z 2 )< 7 w(z T , W(Zi)) = 0 

Since the first equation is non zero, therefore [^] gu(z T ,U(Zi)) f 0. Hence h u (Z 2 ) = 0, which 
means that Vz T , /(Z 2 ; z T ) = 0 i.e. Z 2 ^ 7 ^q( 3) (due to the definition of 7 £q( 3)) which contradicts 
our assumption that Z 2 G 7^©(3) and the proof is complete. ■ 

Remark 1. Zmss A not complete. To show this, we provide a measurable function g(.) / 0 
w’/iere E[g(Z M55 )] = 0/or all z T . Define the centroid (average) of Z MSS by z = z » 

and the center of mass of convex hull of Z MSS by z = j —-—— fc(z M ss) z ^ v w ^ ere C (%mss) 

C ( Z MSS ) 

2 sub-script U is used to differentiate the factorization for sufficient statistics U from the one done for sufficient statistics T 
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represent convex hull of the set, Z MSS , and f zdv is the integral over hyper volume. Since G is 
assumed to be sufficiently large, the distribution of the location of detecting sensors is isotropic 
around the target. Thus, so is the distribution of z and z. Thus E[ z] — E[ z] = z T . In general, 
when dimension of the space is higher than 1, z^z. Consequently, for g (Z MSS ) = z — z we 
have E[g(Z M ss)\ = E[ z — z] = 0 for all z T but z - z ^ 0. 

In one dimensional space, the Zmss simply becomes the maximum and minimum of z, G Z. An 
equivalent problem to this case has been studied in I l32\l and the conclusion is that when R is 
known, Z M ss is not complete. 

Theorem 2. If Z f f], then T(Z) is convex. 

Proof: We draft the proof by mathematical induction. Define Z k = (z 1 ,.., z k ) where k < n. 
Because 2^0, T(Zi) is a ball and hence is convex. If T( Z k ) is convex, so does T{Z k+ f) = 
T(Z k ) fl B R (z k+1 ) because both T(Z k ) and B R (z k+l ) are convex. ■ 


B. Unknown Detection Radius 

In this sub section we assume that R is unknown. This implies that the target’s transmit power 
is not known. Let the polytope C(Z) = A;z, |zj e Z , X/ A* = 1, A, > 0} be the convex hull 

of the locations of the detecting sensors 03ll . Moreover, let W = {w,} c Z be the set of 
locations of comers of C(Z). 

Lemma 2. Any sensors located in C(Z) is a detecting sensor. 


Proof: Since z, e C(Z), it can be represented as z, = X.; A 7 w y where X.y A ; = 1 and 
A j > 0. Thus, 


z i - z T || 


E x 3 w i " z tE A i 

3 j 


]T Xjiwj - z T ) 

3 


< 5Z Aj|I w j- — zxii <J2 x j r 

3 3 


R. 


Theorem 3. W w a sufficient statistic for estimating z T and R. 
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Proof: The probability density function that the detecting sensors will be located at Z is 

f(7 m ( 1 \ w(z) / A(G) - A(B R (z T )) \ N ~ n(Z) n . 

V^(G)J V 4(G) ) ZieB,{(ZT) 


/ 1 \ n(z) /4(G)-4(£ R (0))\ N - n(z) 

U(G)J V 4(G) J Zi 


z T ) 


where n(Z) represents the number of elements of vector Z and O is the origin. Due to Lemma 
[2j if all elements of W are inside B R ( z T ), so is C(Z). Thus all detecting sensors are also inside 
Bf> (zj) and the value of corresponding indicator functions are 1. On the other hand, if any 
Zi G W is not located inside Br{ z t ), / (Z; z T , i?) would be zero regardless of value of indicator 
function for other sensors. Hence, equation Q becomes 


/ (Z; z T , -R) 


i y (z) 

.A(G) J V ^4(G) 


JV-n(Z) 

n ^ z iss fl (z T ) 

ZiSW 


Fisher-Neyman factorization of (|4]) is 


/(Z; z T , f?) = h(Z)g Zl}R (W), 


(4) 


where 

{ fi'z-r.R (VV) = rU eW lz ;eBfl (zT), 

MZ) = (^)’ ,,Z, ( - 4|G, -,1g" (0,, ) A '"’ ,,Z) - 

This implies that, W is a sufficient statistic for estimation of z T and R due to Fisher-Neyman 
factorization theorem Il34l . ■ 

It’s worth noting that for some large R, all corners of W contribute to the possible target 
region. We could not come up with a formal proof for this in l dimensional space but a conjecture 
in two dimensional space is that if it is not, the non-contributing corner would be located in line 
with the other two adjacent corners which contradicts that it is a comer of convex hull. Thus, 
with the same reasoning as proof of theorem [T] W would be a minimal sufficient statistics. Note 
that in this case R, z T are both parameters of estimation. 


Remark 2. Although W is a sufficient statistic for estimation of target location, it is not complete. 
To see this we will provide a measurable function g(.) f 0 where -E[g(W)] = 0 for all z T . 
Denote the centroid (average) ofW by z(W) = ^ Y.,jgw w j> an d the center of mass of C(W') 
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by z (W) = ,■ — 1 - 1 f cm z du. As in RemarkUl it is straight forward to show that E[z] — E[. 

Jc(VV) v 


U 


z = 


zj. On the other hand, in general z/z (except in a one dimensional space). Consequently, for 
g(W) = z(W) — z(W) we have -Z?[g(W)] = E[z — z] = 0/or o// z T but g(W) / 0. 

The one dimensional space is an exception in this regard and has been studied in 021/ . The 
conclusion is that when R is unknown, Z M ss complete in this case. 


Remark 3. Recall that in Theorem [i] we showed that W is a sufficient statistic when R is 
unknown. Hence, this set is also a sufficient statistic for the case when R is known. That is, if 
we define 

Sss = {i\zi e W} 

Zss = 

we have Z MSS C Z ss . Thus, considering that finding the convex hull of detecting sensors is a 
straight forward procedure, Sss can be used as a starting point for building Smss for the case 
when R is known. 


Smss — {k\keSss,T k (Zss)fT(Z S s)}. 


IV. On the Existence of MVU Estimator 

It is easy to show that the likelihood function is discontinuous, and thus, the regularity 
conditions required to obtain Cramer-Rao bound (CRB) does not hold. On the other hand, since 
the minimum sufficient statistics is not complete, Lehmann-Scheffe theorem ll34l cannot be used 
to find the MVU estimator. In this chapter we explore strategies to find the MVU estimator 
directly from its definition while CRB can not be established. The approach lead to interesting 
results regarding the existence of an MVU estimator in the cases when R is known or when it 
is not known. 


A. Known Detection Radius 

Lemma 3. Assuming a target location z T , an obserx’ation Z is a possible observation i.e. 
/(Z; z f) > 0 if and only if z T G T(Z) 
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Proof: The forward direction has already been shown during definition of T. The reverse 
is true because 

z T € T(Z) =>- z T G P| £>p, Zi 

ies 

=>\/i G S, || z T — Zi || < R 

Thus that observation is a possible observation. ■ 

Let us define near points, J\f(T (Z mss )) as all the points whose maximum distance from T (Z mss ) 
are less than R, i.e., 


A/"(T(Z mss )) = {z|z G G & max ||z — zp|| < f?} 

z h &T(Z 

mss) 

in other words if a detecting sensor located inside A^(T(Z mss )), it does not contribute to T and 
is not a member of minimal sufficient statistics. 

Figure |TJa) shows an example of such region. The detecting sensors are demonstrated by small 
empty circles. The target is shown by an empty triangle. The white area is T( Z mss ) region, dark 
gray area illustrates the boundaries of J\f (T(Z mss )). The black circle identifies the detection disk. 
Note that if T(Z mss ) 7^ 0 , Af(T(Z mss )) is not empty because at least Z mss c A/"(T(Z mss )). 

Lemma 4. Af(T(Z mss )) c B R (z T ). 

Proof: We know that z T G T{Z mss ). Thus, according to definition of Af we will have, 

V z P G J\f(T(Z mss )), ||Zp - z r || < R => z P G B r {z t ) 
thus, the proof is complete. ■ 

Theorem 4. Center of mass of T( Z ) is the minimum variance estimator conditioned that at 
least one sensor is detecting and the target is located well inside the sensor deployment region 
such that B‘ 2 r(z t ) C G. 

Proof: Assume that Z mss is the minimal sufficient statistics of Z and Z rass is the vector 
format of that. Then according to lemma [lj T(Z) = T(Z M ss)- 

Let’s assume that an MVU estimator exist for estimation of z T and is represented by 6{ Z). 
Since Z mss is the minimal sufficient statistics, h(Z mss ) = E[6( Z) | Z mss ] should also be an MVU 
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(a) 


(b) 


Fig. 1: T and A f(T) shifts along with Z mss 



Fig. 2: T and A f(T) rotates along with Z mss 


estimator according to Rao-Blackwell theorem lf34l . Moreover, because Z mss can be considered 
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as a function of Z mss , we can restrict our search for MVU estimator to functions of Z mss . In 
other words, if MVU estimator exists we can express it as a function of Z mss i.e. 

9 (^mss) (6) 


On the other hand, the conditions £> 2 r(z t ) C G and at least one sensor would be detecting 
guarantees that for any possible T = T(Z mss ), if Z mss elements shift by some vector d, |d| < R 
then A/”(T(Z mss )) will not reach to the borders of G and its area remains fix. We will calculate 
the variance of location estimation by conditioning over the number of detecting sensors, n, 
and over the number of sensors forming the minimal sufficient statistics denoted by k. That’s 
because the incident of each ( n,k ) partitions the random space to disjoint subspace. 

Let us represents (zi,... z k ) by Z k . Considering that below expectation is symmetrical regardless 
of the order of the sensors forming the minimal sufficient statistics, we assume that the sensors 
belong to minimum sufficient statistics have indexes 1,.., k and generalize the result by including 
a binomial coefficient. The assumption that the minimal sufficient statistics are located at index 
1 dictates that other detecting sensors should be located in J\f(T(Z k )). Therefore, the 
expectation of square error of the estimation conditioned that at least one sensor is detecting can 
be calculated as following: 


771 I ~ 12] 

b[\zj — zj\ 


/ .4(G) - A(B R (0)) \ N ~ n (N\ " (n 

k\ A ( G ) ) w k\ k . 


>B r ( zt) Jb r (z t) ^AT(T(Z k )) JJ\f(T(Z-k)) 


k 

1 \ n 


n—k 
12 


.4(G)/ 1/ ( z k) - z t|" !z k =MSS(Z k ) 
dz n ..dz k - k idz k . .dz\ 


where N represents the total number of sensors, L dz 



dxi..dxi is the integral over 


i 

the multi dimensional volume, Q, and j is the probability that N — n sensors 

located outside the detecting disk. The indicator function lz k =MSS(z k ) make the integrand zero 
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whenever Z k is not the minimal sufficient statistics. 


N 


£[Kt-z t ,-i - E[ n JE w 


N\ ” (n\ (4(G) - A(S s (0))) 


N—n 


.4(G) 


N 


!Br{ zt) JBr( zt) W(Z k )) JM(T(Z k )) 


dz„..dz 




n—k 

|2 


lz k =MSS(Z k ) Is (Zk) — Z T </Zfc..dZi 

E f^) E ^ ( ' A(G) " • 4(Ss( ° )))W ~“ 


n=l 


n 


w 


.4(G) 


AT 


^h( z t) ^r( z t) 


(A (Af(T(Z k )))) n - k \g (Z k ) - z T | 2 l Zk =MSS(Z k ) 
d,z k ..dzi 

Now, if we change the variables according to the following rules 

tj = 0 

t 2 = z 2 — Zi 


f k ^fc ^1 


and assuming T k = (0,t 2 ..,t fc ) and zil k = (zi, ...zi), then (zi, ..,z fc ) = zil k + T k . Moreover, 

k 

it is reasonable to assume that the MVU estimator, g, shifts in space whenever the the input 
data shifts i.e. g (zil k + T k ) = g (T k ) + zi cause otherwise the g would be dependent on the 
selection of the origin. 


771 M ~ I 2' 

h/ [|z j — Zt | 


f ( N ) V ( n ) (AG) - A(B R (Q))) N - 

h\n)kV) AG) N 


j S(zi) 


(zi) 


(A (A/”(T( Zl l k + T k )))) 


n—k ■ 


*-z T GT(zil k +T k ) 


X 


/c-1 


I# (T k ) + Zi — Z T |" l(zil k +T k )=JW55(zil k +T k )<ytfc-- ( yt2(iZi 


where the interval of integrals for elements of T k are the shifted versions of Br( z t ) by z± i.e. 
S ( z ^) = BidZj) — z 1 . Moreover, the indicator function, l Zxe r ( z 1 i k + T k )^ make the integrand zero 
whenever T k is an impossible (invalid) combination. 

The next equation uses the fact that when £> 2 _r(z t ) C G, shifting T k by zi will not change 
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the area of A{N(T)) i.e. A(J\f(T( T k + zylij)) = *4.(Af(T(T k ))) as illustrated in figures 
[TJa) and [IJb). That’s because according to lemma (|4|), A/"(T(T k + zil k )) C Br( z t ). It is 
also worth noting that the result will not change if we substitute the inner integral intervals, 
§(z^), with B 2R (0) = Br(O) © B r (0) (where © represents the Minkowski sum 0510 because 
§(z | j c B 2 r(0) and the indicator functions, 1, guarantees to make the integrands zero whenever 
tj ^ §(zi) & ti G B 2 r(0 ) for any i £ (2, k}. Thus, we will have 


N 


N' 


e[ |z T - z T | ] = E L E 


'n\ (A{G)-A(B r ( 0))f~ n 


=i W t w 

r / ../ 

B r { z t ) Jb 2R { O) Jb 2R {O) 


A(G) N 

(A(AT(T( T k )))) n ~ k 


-z T eT(zii k +T k ) 


X 


-V- 

k-1 


|g (T k ) + Zi — Z T | 2 l(zil k +T k )=MSS(zil k +T k )^tfc--^t2^Zi 

On the other hand, it is obvious from the mechanism of building T that T(zil k + T k ) would 
be the shifted version of T(T k ) by zi, and similarly M,S',S'(zi l k + T k ) would be the shifted 
version of MSS( T k ) by zi. Thus, 


N 


e[\z t - zri ] = El L EZ 


72—1 


N\ " kn\ (Aj G) - .4(^(0))) 
n ) ^ \kj A(G) N 
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" -V-' 

fc-1 

| <7 (T k ) + Zi — Z T | 2 lT k =ArSS(T k )^tfc--^t2^Zi 


Now with a change of variable t' = z T - z 1; we have 
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N > 


E [ | Zj - z T | ] = E L E 


'n\ (^(G)-^(S r (0))) 
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n= 1 


n 


t w 
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7zt—B_r(zt) JB 2 r(Q) JB 2 r( O) 


fc-1 


I# (T k ) — t/| 2 lt' e r(T k )rftfc••rft2c?t , 


Now because the intervals of the inner integrals are independent of t', we can change the 
order of outer integral with the inner ones. Moreover, because the first integral interval is taken 
over the hyper volume, we know that I Zt _ Br{zt) - dt ' = / Bij(ZT) _ ZT -d(-t') = J Br{0) ..d(-t') = 
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jj 3 lt (o) ■■d(t '), thus we can reorganize the integral as: 
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E [ I Z T - Zx| ] = E L E 


N\ " / n\ (A(G)-A(B R m) 
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n= 1 
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JV 


/ •• / (^(•^(7 _ (Tk)))) n fc lT k =MSS(T k ) X 

JB 2 r(0) JB 2 r{0) 

-v-- 

fc -1 


/ lt'er(T k ) |p (Tk) — t 7 | fit' j dtk-.dt-i 
>Br( O) ) 

because at least one sensor is assumed to be detecting T(T k ) C B R (0). Therefore, the last 
integral can be written over T(T k ) i.e., 

i?[|zt— z T i 2 ] = £ pi ± h) ^ ~ 


N-n 


n= 1 


n 


\ k J 


A(G) 


N 


[ ■■ / (-4(A/"(T(T k )))) n k l Tk =Mss(T k ) x 

Jb 2R {o) Jb 2R { o) 


fc-i 


'r(T k ) 


|p (T k ) - t'| fit') dt fe ..dt 


(7) 


which is indicating that to make the expectation minimum, g (T k ) should be the center of mass 
of T(T k ), i.e. 


9 (T k ) = CM (T(T k )) 


( 8 ) 


where CM denotes the center of mass. Now that the answer to the minimization problem is 
known, it is easy to verify that g (T k ) = CM(T( T k )) makes the variance minimum because 
g (T k ) is independent from t' and the only point in space which has minimum average distance 
from all t' <E T(T k ) is CM (T(T k )). In other words, it makes the inner integral minimum and 
g (T k ) does not show up in any other parts of the expression. It’s worth mentioning here that 
although center of mass of a possible target region has been used as a locali z er in literature 
ll36l . Il371l . flMS, it appears that it has never been proved to be the MVU estimator. 

Furthermore, similar to previous derivation we can show that g (T k ) = CM (T(T k )) is unbiased 
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as following: 

E[ Zt - z T ] 


f ^ V H ( ^ (G) - W»))^ 

.4(g)" 

/ •• / (-4(A/"(T(T k ))))" fc l Tk =Mss(T k ) x 

•AiWO) JB2r(0) 

S --V- 

fc -1 

f / lt'eT(T k ) (p (T k ) — t 7 ) fitfc..fit 2 
V'iT(T k ) ) 

f ^ V M (4(G)-^OMO)))”-" 

iwiW 4(G)- 7 

/ •• [ (4(A/"(T(T k ))))" fc lT k =Mss(T k ) x 

JB2r{0) JB2r{ O) 
fc -1 

(0) dt fc ..rft 2 
0 


(9) 


Thus, MVU estimator exists and is unique. From ([8]) and shifting sensor locations by z ls we 
conclude that for any Z e 7^e(3), 


g (Z) = CM (T(Z)) 


( 10 ) 


is the minimum variance unbiased estimator for z T . 


/>’. Unknown Detection Radius 

For this problem the estimation parameters are [z T , i?| and we know that Z ss is a sufficient 
statistics. We consider two detecting range Ri and R 2 . Let us define T(Z,R) = f) ieS Bn(zi). 
Based on previous section analysis, g\ (Z) = CM (T(Z, Ri)) is the unique unbiased minimum 
variance estimator for z T when R = R x and g 2 (T) = CM ( mathcalT(Z, R 2 )) is the unique 
unbiased minimum variance estimator for z T when R = R 2 . Note that no matter what the R 
parameter is, g 2 and g x remains unbiased for location estimation because the probability density 
function of observation is isotropic and if the observation rotates, so do g 2 and g \. These two 
functions can be different for a specific realization of Z as demonstrated in an example in figure 
[3j Therefore, the MVU estimator does not exist for this case because for different parameter, R, 
the minimum variance unbiased estimators are different. 
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Fig. 3: Illustration of T(Z, R ) and corresponding CM for R = 1.05 and R = 1.6 


V. Sub-Optimum Estimators in Two Dimensional Space 


A. When Detecting Radius is Known 


Although we proved in theorem [2] that T(Z) is a convex shape, it is not an easy shape 
to work with. Not only its visualization but also finding the center of mass of that shape is 
computationally costly. An intuitive approach to decrease the complexity is to find the center of 
mass of the convex hull of the comers of T(Z) as illustrated in figure [4j 
The center of mass of this shape can be found easily in two dimensional space by triangulating 
it through a triangulation algorithm such as the ones discovered by Euler, Fournier or Toussaint 
ma-ffl. Assume the result of triangulation to be q triangles each represented by its corners 
Ai,Bi,Ci as A a . )B . a , then 


CM{SCT{ Z)) = 


A-i + Bi + Ci 


X 


( 11 ) 


would be a sub-optimum estimator for z T where SCT{ Z) is the convex polygon generated by 
connecting the consecutive corners of T(Z) by straight lines, and A+A+c, j s p lc ccn t cr of mass 
for the ith triangle 
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Fig. 4: An illustration of SCT{ Z) 


B. When Detecting Radius is Unknown 

As we discussed in previous section, the MVU estimator does not exist for this case. Still 
the MVU of the case when R is known act as a Clairvoyant estimator and its performance 
act as the lower bound for all estimators who lacks the knowledge of R. We may notice that 
all estimators provided in |[20l can be employed when R is unknown. In the next section we 
will see the simulation results and we will find out that center of MEC performance follows the 
Clairvoyant estimator closely. So it may be used as a sub-optimum estimator in this case. 


VI. Simulation Results 

A set of simulations are performed in two dimensional space to verify that for known R, 
the center of mass of T(Z) is the minimum variance estimator among other famous estimators. 


’A Clairvoyant estimator is referred to an estimator with advanced knowledge of some parameters of estimation as if it is 
provided by a genie. 
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We assume that in each trial N sensors are dispersed randomly in a rectangular region A = 
100m x 100m, centered at the origin, where the target is located. Number of sensors depend 
on the density of sensor deployment, p and is assumed to be N = \_p x A\. Any sensors 
located within the R distance of the origin would be a detecting sensor and vice versa. In f[20ll . a 
number of heuristic algorithm introduced for estimation of target location when only the location 
of detecting sensors will be reported. In these methods, the target location will be estimated as 
the Steiner center, the center of the minimum enclosing rectangle or the center of the minimum 
enclosing circle for the locations of detecting sensors denoted by Steiner center, center of MER, 
and center of MEC respectively. For the purpose of comparison, the center of mass of possible 
target region based on the detecting sensors, CM (T( Z)), have been considered along with the 
above mentioned heuristic methods. An extension of a grid base algorithm has been used to find 
the center of mass of possible target region. In addition, Delaunay algorithm readily available 
in Matlab is employed to implement the triangulation used in sub-optimum estimator PTl . [[421. 
Ignoring the trials that result in no detecting sensor, mean square error of these methods have 
been calculated for each case. 

Figures [5] and [6] depict Mean Square Error (MSE) of the above mentioned estimators versus 
density when R is fixed and is equal to 1 and 5 respectively. The results have been averaged over 
8000 trials. As can be seen in the graphs, the center of mass of T(Z) beats heuristic estimators 
as density increases. When the density is small, the likelihood of trials with only one or two 
sensors detecting is high which in these cases, all methods have identical estimates. It is also 
noticeable that for large densities, the sub-optimum estimator follows the CM (T(Z)) closely. 

Figure [7] and [8] depict MSE of the above mentioned estimators versus R when density of 
sensor deployment, p, is fixed and is equal to 1 and 4 respectively. The results have been 
averaged over 2000 trials. As is clear from these graphs, the center of mass of T{ Z) beats the 
heuristic estimators again as R increases. 

VII. Conclusion 

This paper investigated the existence of an MVU location estimator when noise free detectors 
are deployed around a target and only the detecting sensors report their locations to a fusion 
center. It is proven mathematically that when the radius of detection is known and the target 
is located well inside the sensor deployment region, the center of mass of the possible target 
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Fig. 5: MSE versus density for R=1 


region is the MVU estimator and when the radius of detection is not known the MVU estimator 
does not exist. Moreover, the minimal sufficient statistics of the detecting sensors are derived 
both when the radius of detection is known and when it is not known. In addition, a set 
of simulations is performed to compare the performance of the MVU estimator with various 
heuristic estimators. Finally, a sub-optimum estimators introduced, which is computationally 
less complex and the processing cost is independent of its resolution. It is also shown that when 
the density of deployment sensor is increased the performance of the sub-optimum estimator 
approaches the MVU performance. 
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